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This Letter describes measurements of inclusive W(— > ev) + n jet cross sections (n =1-4), pre- 
sented as total inclusive cross sections and differentially in the n th jet transverse momentum. The 
measurements are made using data corresponding to an integrated luminosity of 4.2 fb _1 collected 
by the DO detector at the Fermilab Tevatron Collider, and achieve considerably smaller uncertainties 
on W+jets production cross sections than previous measurements. The measurements are compared 
to next-to-leading order perturbative QCD (pQCD) calculations in the n =1-3 jet multiplicity bins 
and to leading order pQCD calculations in the 4-jet bin. The measurements are generally in agree- 
ment with pQCD calculations, although certain regions of phase space are identified where these 
predictions could better match the data. 

PACS numbers: 12.38.Bx, 13.85. Qk, 14.70.Fm 



Measurements of vector boson plus jet production are 
fundamental tests of perturbative quantum chromody- 
namics (pQCD), the theory describing the strong inter- 
action. In addition to providing a test of pQCD at high 
momentum scales, VF+jets production can be the dom- 
inant background in measurements of single top quark 
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and tt production as well as in searches for the stan- 
dard model Higgs boson and for physics beyond the stan- 
dard model. Theoretical uncertainties on the production 
rates and kinematics introduce limitations in our ability 
to identify new physics signals. Therefore, it is crucial 
to make precision measurements of VF+jets production 
at the Fermilab Tevatron Collider and the CERN Large 
Hadron Collider in order to constrain these backgrounds. 
We present new measurements of VF+jets cross sections 
with a data sample more than ten times larger than that 
used in previous measurements [l[ , allowing the first de- 
tailed study of W + 4 jet production. The previous mea- 
surements have been used extensively in testing and tun- 
ing theoretical models of W boson production 0-0] ■ 
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The strategy employed for this measurement is based 
on those used in the DO Z+jet cross section [B| and Z bo- 
son p T [6] publications. We select a high purity sample of 
Vy+jets events and the results are corrected to the "par- 
ticle level," which includes energy from stable particles, 
the underlying event, muons, and neutrinos, as defined in 
Ref. 0- This procedure corrects a measured observable 
back to the particle level observable, correcting for the ef- 
fect of finite experimental resolution, detector response, 
acceptance, and efficiencies. 

These measurements use a sample of W(— > ev) + n 
jet candidate events corresponding to an integrated lu- 
minosity of 4.2 fb _1 collected with the DO detector in 
Run II of the Fermilab Tevatron Collider. The DO de- 
tector consists of a central tracking system, comprising a 
silicon microstrip tracker and a fiber tracker, both within 
an approximately 2 T axial magnetic field. These com- 
ponents are used primarily to identify the location of the 
pp interaction vertex and the electron produced in the 
decay of the W boson candidate. Outside of the track- 
ing system, a liquid-argon and uranium calorimeter is 
divided into a central section and two end sections that 
are used to identify electromagnetic and hadronic show- 
ers. A detailed description of the DO detector can be 
found in Ref. Q. 

The data were collected using a suite of electron and 
electron+jet triggers. The lowest electron transverse en- 
ergy threshold in the electron suite is 22 GeV, and the 
electron threshold for the e+jets triggers is 15 GeV. The 
combination of the triggers used provides > 97% trig- 
ger efficiency for electrons with transverse energy above 
26 GeV. The efficiency in the turn on region below this 
energy threshold is evaluated using unbiased data sam- 
ples and a corresponding scale factor is then applied to 
the MC simulation. 

The events were then processed through the DO re- 
construction program which identifies jet and W boson 
candidates. Jets are identified with the DO midpoint cone 
algorithm [9j, which uses a cone of radius 1Z = 0.5 (dis- 
tance in 77— 4> space 10]) to cluster calorimeter cells. The 
electromagnetic fraction of the jet energy is required to 
be below 0.95 to reject electrons and above 0.05 to sup- 
press jets dominated by noise. Jets with a large fraction 
of their energy deposited in the coarse hadronic layers 
of the calorimeter are also rejected due to noise typi- 
cal in those layers. To minimize background from jet 
candidates arising from noise in the precision readout of 
the calorimeter, confirmation from the readout system of 
the first level trigger is required for reconstructed jets. 
Jets matched to loose electrons with pr > 20 GeV and 
Ai?(e,jet) < 0.5 are also rejected. Jets are corrected 
for calorimeter response, instrumental and out-of-cone 
showering effects, and additional energy deposits in the 
calorimeter that arise from detector noise and pile-up 
from multiple interactions and different beam crossings. 
These jet energy scale corrections [llj are determined us- 



ing transverse momentum imbalance in 7 + jet events, 
where the electromagnetic calorimeter response is cali- 
brated using Z/j* — > e + e~ events. Jets are required 
to have at least two tracks that point to their associ- 
ated pp vertex. Energies of jets containing muons are 
corrected with the measured muon momentum after ac- 
counting for the typical energy deposited by a minimum 
ionizing particle. Jets are ordered in decreasing trans- 
verse momentum and we call the jet with the highest 
transverse momentum "leading." Electrons are identified 
as clusters of calorimeter cells in which 95% of the energy 
in the shower is deposited in the electromagnetic (EM) 
section. The electron candidates must be isolated from 
other calorimeter energy deposits, have spatial distribu- 
tions consistent with those expected for electron show- 
ers, and the event must contain a reconstructed track 
matched to the EM shower that is isolated from other 
tracks. Isolation from energy deposited by hadrons is im- 
posed by requiring (E tot - E cm )/E cm < 0.15, where E tot 
(E cm ) is the total (electromagnetic) energy in a cone of 
radius 1Z — 0.4 (1Z — 0.2). Events with a second isolated 
electron (with px > 15 GeV) are removed to suppress 
the background due to Z boson and Drell-Yan produc- 
tion. The missing transverse energy in the event is cal- 
culated as the vector sum of the calorimeter cell energies 
and is corrected for the presence of any muons. Because 
the longitudinal component of the momentum of the neu- 
trino is not measured, the measured properties of the W 
boson candidates are limited to their transverse energy, 
E)p , and transverse mass, defined as 
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where is the magnitude of the missing transverse en- 
ergy vector, pip is the transverse momentum of the elec- 
tron, and p e x and p y and p" y ) are the magnitude of the 
x and y components of the electron's momentum (missing 
transverse energy) respectively. 

The following requirements are used in order to sup- 
press background while maintaining high efficiency for 
events in which a W boson is produced: Pj. > 15 GeV 
and electron pseudorapidity |?7 e | < 1.1, $ T > 20 GeV, 
> 40 GeV, jet transverse momentum p 3 ^, > 20 GeV 
and rapidity |y jct | < 3.2, ATZ = ^{Acf)) 2 + (A77) 2 be- 
tween the electron and the nearest jet > 0.5, and the z 
component of the pp interaction vertex is restricted to 
I z v t x I < 60 cm [l(| ■ Events must have a reconstructed pp 
interaction vertex, containing at least three associated 
tracks. This pp interaction vertex is required to be less 
than 1 cm away in the coordinate along the beam line 
from the extrapolated electron track. 

After these requirements, W^(-l-jets) events dominate 
the data sample but there are backgrounds from Z+jets, 
W{— > tv — > ew)+jets, tt, diboson, single top quarks, 
and multijet events. We simulate the W/Z+jets and 

processes with ALPGEN [12( interfaced with PYTHIA [13 1 
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for the simulation of initial and final state radiation and 
for parton hadronization. The pythia generator is used 
to simulate diboson production, while production of sin- 
gle top quarks is simulated with the COMPHEP 14| gen- 
erator interfaced with PYTHIA. The cross sections for 
W/Z+jet production are taken from ALPGEN, corrected 
with a constant multiplicative factor to match the in- 
clusive W/Z+jet cross sections calculated at NLO fljj . 
Additional corrections are applied to events containing 
W/Z bosons plus heavy flavor jets, to match the predic- 
tions of NLO QCD calculations. Events from randomly 
chosen beam crossings, with the same instantaneous lu- 
minosity profile as the data, are overlaid on the simu- 
lated events to reproduce the effect of multiple pp inter- 
actions and detector noise. All simulated samples are 
passed through the DO detector simulation and then re- 
constructed in the same way as the data. The estimated 
fraction of the data sample that is due to processes other 
than V^+jets ranges within 2-40%. Leptonic background 
from W(— > tv — > ez/i/)+jets processes represents approx- 
imately 5-8% of all reconstructed MK+jets events, and 
the fraction of background due to top quark production 
ranges within to 7% (16%) in the one (two) jet multi- 
plicity bin, 5-40% in the three jet bin and 20-60% in the 
four jet bin (with the extremes only being reached at the 
highest jet px bins in all cases). 

In multijet events, there is a small but non-negligible 
chance that a jet may be misidentified as an electron and 
then the event may pass all selection criteria. As the 
multijet cross section is large, the contribution from such 
instances of fake-electron events to the measured distri- 
butions must be taken into account. To determine the 
number and kinematic distributions of such events, we 
use the data-driven method described in Ref. .16] because 
the estimation of this background from Monte Carlo sim- 
ulations is not reliable. This approach uses data in a 
control region that has no overlap with the signal selec- 
tion to determine the differential distribution and overall 
normalization of the multijet distributions. 

The total background contribution is subtracted from 
the data in each bin of the p£ distribution. After back- 
ground subtraction, the data are corrected for detector 
resolution effects using a regularized inversion of the reso- 
lution matrix as implemented in the program GURU [l7| , 
with ensemble testing used to derive statistical uncer- 
tainties and unfolding biases. This method is described 
in detail in Ref. [6j. We have chosen the matrix unfold- 
ing approach over the traditional bin-by-bin correction 
method because of non-negligible bin migration effects 
in the pt£ variable and because the matrix unfolding 
method provides improved estimation of the uncertain- 
ties of the measurement. 

To evaluate statistical uncertainties on the unfolded 
distributions, as well as systematic biases and uncertain- 
ties, we build ensembles using ALPGEN+PYTHIA signal 
events that have the same statistical fluctuations as the 



data sample. The ensembles are reweighted to accurately 
describe the kinematics of the unfolded jet px- Five hun- 
dred ensembles are created and unfolded in the same 
manner as the data and are in-turn compared to their 
corresponding generator-level distributions. The resid- 
ual differences between the generator-level and unfolded 
measurement in each bin, for each ensemble, are deter- 
mined and fitted with a Gaussian function. The mean 
offset of the distribution is used to construct an unfold- 
ing bias correction to be applied to the data, while the 
larger of the root mean square and the Gaussian width 
is assigned as the statistical uncertainty associated with 
that bin in the unfolded distribution. The unfolding bias 
correction is small, generally 0.5-2%, and always much 
smaller than the statistical uncertainty in the bin. Over- 
all, the statistical uncertainties are within 1-17%, de- 
pending on jet multiplicity and jet pr bin. 

The systematic uncertainties affecting this measure- 
ment can be divided into three types: those related to 
the knowledge of the detector response, those related to 
the background modeling and those associated with the 
unfolding method itself The systematic uncertainties re- 
lated to the modeling of the detector response and their 
effect on the final cross sections arise from the calibration 
of the jet energy scale [3-16%], from the measurements of 
the jet energy resolution [0.1-17%], the jet identification 
efficiency [0.3-4%], the jet-track matching requirement 
[1-11%], the trigger efficiency [1-4%], the electron iden- 
tification efficiency [4-5%], and the uncertainty in the 
luminosity determination [6.1%]. We determine the sys- 
tematic uncertainty for all these sources apart from the 
latter two using the ALPGEN+PYTHIA ensembles. The 
relevant variables in all events are varied within their sys- 
tematic uncertainties, resulting in new signal templates 
and new migration matrices. The nominal ensembles 
(which look and behave as our reconstructed data dis- 
tributions) are again unfolded but this time with inputs 
to GURU replaced with the systematic-shifted samples 
As expected, it is found that the statistical uncertain- 
ties from the shifted residual distributions are largely in- 
sensitive to changes in the detector response, but the 
unfolding bias can vary significantly. The change in the 
bias from the nominal to shifted ensembles is attributed 
to the systematic uncertainty in the unfolded data distri- 
butions. All differential cross section measurements are 
normalized to the measured inclusive W boson cross sec- 
tion, resulting in a complete (partial) cancellation of the 
systematic uncertainties due to luminosity (trigger and 
electron identification efficiencies). The dominant uncer- 
tainties due to jet energy scale and jet energy resolu- 
tion are correlated bin-to- bin (and between jet spectra), 
the uncertainties due to the jet-track matching require- 
ment and electron identification efficiency are partially 
correlated. All other uncertainties are considered to be 
uncorrelated. The correlation of systematic uncertain- 
ties between jet multiplicity bins are taken into account 
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when normalizing the differential cross section spectra 
and in determining the uncertainties on measurement of 
the (T n /cr n _i inclusive cross section ratios. 

The remaining sources of systematic uncertainty are 
the normalization and differential distributions of the 
multijet background [0.1-4%], the uncertainty due to the 
electron final state radiation at particle level (<1%), un- 
certainties associated with the unfolding method (<1%) 
and the theoretical uncertainty on the tt cross section. In 
some regions of phase space (at high pr in the three and 
four jet multiplicity bins) the data sample is dominated 
by tt production. In these regions the ~ 8% uncertainty 
in the tt cross section translates into an uncertainty of 
up to 19% in the tt subtracted VT+jets signal. Uncer- 
tainties due to the unfolding procedure come from the 
uncertainty on the derivation of the unfolding bias used 
to correct the unfolded spectra, and from the change of 
the final result when this is obtained repeating the un- 
folding procedure with a data-derived reweighting of the 
MC inputs to GURU in order to account for mismodeling 
effects present in the Monte Carlo predictions. 

As in the case of the differential cross section measure- 
ments, the inclusive W(— ¥ e^)+jets production cross sec- 
tions are normalized to the measured inclusive W — > ev 
cross section. This normalization reduces (or cancels) 
systematic uncertainties and provides sensitivity to the 
shape of the distribution in comparisons to Monte Carlo 
and theoretical predictions. The events passing the selec- 
tion requirements arc well described by the Monte Carlo 
predictions and the sample is dominated (> 99.8%) by 
the inclusive production of W events. The total inclu- 
sive W boson cross section within the kinematic accep- 
tance is measured to be aw — 1097 ± l(stat) +g 9 (syst) 
± 67(lumi) pb. This number is used to normalize the 
differential cross section results. 

Recent theoretical work 0, [H| nas extended the avail- 
ability of predictions up to W+3 jet events at NLO. 
Although there has also been a recent calculation of 
W+4 jet production at NLO for pp collisions at y/s = 7 
(or 14) TeV these predictions are not available for 
the Tevatron, and comparisons with theory are there- 
fore limited to LO for W+A jet production. In this 
analysis, we use the interfaced BLACKHAT+SHERPA (2f3 | 
and rocket+MCFM [2l|, |22| programs as the main 



sources for theoretical predictions of W+jets production. 
The MCFM calculations employ version 6.0 of the pro- 
gram. BLACKHAT and rocket are parton level genera- 
tors which incorporate NLO QCD calculations with up 
to 3 final state jets. They provide parton level jets corre- 
sponding to the hard partons, but they do not include the 
underlying event or hadronization effects. We compare 
both theory predictions to our measured cross sections, 
in order to determine the differences that arise from theo- 
retical choices made in the calculations, such as the choice 
of renormalization and factorization scales, and in order 
to explore the uncertainties inherent in these predictions. 



The blackhat+SHERPA program employs the renor- 
malization (hr) and factorization (jxf) scale fi = fxp = 
fiR = \H' T , where H' T is the scalar sum of the parton 
and W transverse energies. blackhat+SHERPA docs 
not provide cross sections using the DO midpoint jet algo- 
rithm, but instead uses the SISCONE Q algorithm with 
split-merge parameter / = 0.5 and cone radius 1Z = 0.5. 
In order to keep all the theory predictions on the same 
footing, we therefore show the blackhat+SHERPA and 
ROCKET+MCFM predictions using the SISCONE jet algo- 
rithm. The effect of differences in the theoretical predic- 
tions produced with different jet algorithms was found to 
be approximately one order of magnitude smaller than 
the scale uncertainties in all jet multiplicity bins, and so 
is considered to have negligible impact on the interpreta- 
tion of the theory /data comparison. The choice made by 

the rocket+mcfm authors is fj, = J M w + j(T,pi et ) 2 

(in the 2, 3, and 4- jet bins), summing over the four- 
momenta of all jets in the event, where Mw is the mass of 
the W boson. This scale choice was suggested in Ref. [24[ 
because it sums large logarithms in the calculation to all 
orders. In the 1-jet bin, a slightly modified choice of 

fx = J M w + (p J y*) 2 is used. This is due to the fact that 
in the 1-jet bin, the NLO calculation includes diagrams 
with an extra hard (real) emission or virtual loop correc- 
tions. For the Born and virtual loop diagrams, the only 
hard scale is Mw, due to the single massless jet balancing 
the W boson. However, in the case of diagrams with an 
extra hard emission, the two final state partons can be 
combined into one massive jet by the jet reconstruction 
algorithm increasing the scale of the real contributions, 
which generally increase the cross section. As a result, 
the real diagrams are evaluated with a coupling that is 
smaller, due to the running of a s , than the virtual di- 
agrams, which leads to a prediction of the NLO cross 
section that is too low. Both theory calculations use the 
MSTW2008 parton density function (PDF) 0, where 
the LO (NLO) cross section calculation is matched to the 
LO (NLO) PDF. The uncertainties on the theory predic- 
tions are estimated by multiplying \i by factors of 2 and 
0.5. 

Fixed-order pQCD predictions provide only a parton- 
level prediction which is not immediately comparable 
to the unfolded data. Additional corrections must be 
applied to propagate the fixed order predictions to the 
particle level. The two effects which contribute to 
this parton-to-particle correction are hadronization of 
the final state partons and the presence of the un- 
derlying event. These corrections (referred to collec- 
tively as hadronization corrections) are obtained with the 
SHERPA MC program [3], which employs the CTEQ6.6 
PDF set [2^. The corrections are generally around 10%, 
but can be as large as 25% at high pj?. The parton 
level cross sections are computed with the SISCONE jet 
finding algorithm, while the particle level predictions are 
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FIG. 1: (a) Total inclusive n-jet cross sections cr n = a(W( 
ev) + > n jet 



pf > 20 GeV) as a function of inclu- 



sive jet multiplicity, (b) the ratio of the theory predic- 
tions to the measurements, and (c) a n /a n - 1 ratios for data, 
blackhat+sherpa and rocket+mcfm. Error bars on data 
points represent combined statistical and systematic uncer- 
tainties on measured cross sections. The uncertainties on the 
theory points in (a) and (c) and the hashed areas in (b) rep- 
resent the theoretical uncertainty arising from the choice of 
renormalization and factorization scale. In (b) the error bars 
on the points represent the data uncertainties. 



computed with the DO midpoint cone algorithm, in or- 
der to account for the difference in jet algorithm between 
the data and the pQCD predictions. The impact of fold- 
ing the correction for the jet algorithm into the overall 
hadronization correction is small, and approximately an 
order of magnitude smaller than the theoretical scale un- 
certainties in size. All inclusive and differential pQCD 
predictions have the hadronization corrections applied to 
them. We p rovide the tables of the hadronization cor- 
rections [29j so that future pQCD calculations can be 
compared to the data on the same terms. The quoted 
uncertainty on these corrections is purely statistical. 

Fig. QJa) shows the absolute inclusive W + n jet 
cross sections for each jet multiplicity considered, com- 
pared with the LO and NLO theoretical predictions from 
BLACKHAT+SHERPA and rocket+mcfm, where both 
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FIG. 2: Measured W+n jet differential cross section as a func- 
tion of n th jet pr for jet multiplicities n =1-4, normalized to 
the inclusive W — > ev cross section. W+l jet inclusive spec- 
tra are shown by the top curve, the W+4 jet inclusive spectra 
by the bottom curve. The measurements are compared to the 
fixed-order NLO predictions for the jet multiplicities n =1-3, 
and to LO predictions for n = 4. 



are corrected for hadronization effects. Fig. [ljb) shows 
the ratio of theory to data. Good agreement is observed 
between data and the NLO theory predictions, except for 
the 1-jet bin, where the NLO prediction presents a slight 
excess with respect to the data. Fig. [IJc) shows the mea- 
surement of the cr n /(7 n -i inclusive cross section ratio as 
a function of inclusive jet multiplicity for n=l-4 in com- 
parison to predictions of this ratio from LO and NLO 
calculations. Here, the theoretical uncertainty takes the 
correlations of the scale choice between the n and n — 1 
multiplicity bins into account. The data uncertainties 
are also calculated from the relative uncertainties on the 
two cross sections, but with partial or total cancellation 
of systematic uncertainties due to electron identification, 
trigger, and luminosity. The uncertainties due to the 
jet corrections are correlated between bins and are ac- 
counted for. The total uncertainties on the measurement 
presented throughout this paper are comparable to the 
scale uncertainties on the predictions at NLO. Tables of 
the measured and theoretical cross sections and their un- 
certainties are given in the appendix to this paper. 

The unfolded differential data cross sections (multi- 
plied by the branching fraction of the W — > ev decay) 
for each jet multiplicity are shown in Fig. [21 The data 
are normalized by the measured inclusive W boson cross 
section in all jet multiplicity bins, which reduces the un- 
certainties in the measurement because of cancellation 
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of some systematic uncertainties. The data spectra are 
compared to the predictions from ROCKET+mcfm and 
BLACKHAT+SHERPA (again normalized by their respec- 
tive inclusive W boson cross sections and corrected for 
hadronization effects). The theory is able to describe the 
data throughout the pj, spectra for all multiplicities, al- 
though a detailed comparison is best made by examining 
the ratios of theory to data. Each data point is placed 
at the pt value where the theoretical differential cross 
section is equal to the average cross section within the 
bin 



27|. 



The ratio of the theory predictions to the unfolded dif- 
ferential data cross sections are shown in Fig. [3] Each of 
the data and theory cross sections is normalized to its re- 
spective inclusive W boson production cross section. In 
the inclusive W+l jet bin [Fig.[3ja)], the data uncertain- 
ties vary by 4-14%, but for most jet transverse momenta 
these uncertainties are smaller than the theoretical uncer- 
tainties. The data agree well with both NLO theory cal- 
culations, although the theoretical prediction is slightly 
higher than the data at low pji . The inclusive W+2 jet 
bin results are shown in Fig. EJb). The measured un- 
certainties vary by 5-20% and are similar to those of the 
1-jet bin. The blackhat+sherpa and rocket+mcfm 
predictions are in good agreement with the data every- 
where. In Fig. EJc), the ratio of W+3 jet pQCD predic- 
tions to the differential cross sections are shown. The re- 
sults of NLO predictions are below the data at high p^ , 
but still consistent within uncertainties. In Fig. [3][d), 
the differential cross section measurement of W+4 jets is 
shown as a ratio to the LO pQCD prediction. The theory 
prediction can reproduce the data, albeit with large un- 
certainties. Theoretical cross-sections at LO suffer from 
strong dependence on the choice of renormalization and 
factorization scales, in part due to large logarithmic cor- 
rections and higher order contributions. The significant 
reduction of the scale uncertainty at NLO compared to 
the same uncertainty at LO is an indication that the size 
of the NNLO corrections is small. An NLO prediction 
for this final state is necessary to make a more robust 
comparison. 

In summary, W+n jet inclusive cross sections for n — 
1,2,3 and 4 jets have been measured using 4.2 fb _1 of 
integrated luminosity collected by the DO detector. The 
measurements include the total inclusive cross section for 
each jet multiplicity and differential cross sections as a 
function of the n th jet pt- These measurements rep- 
resent a test of pQCD com plem entary to the extensive 
DO Z+jets measurements [5I l28l l30l| . The measured cross 
sections improve on the measurement by CDF 1] by in- 
cluding W+4 jet differential cross sections, by substan- 
tially improving the uncertainties on differential cross 
sections in all jet multiplicities, and by performing the 
first comparison with NLO W+3 jet cross section predic- 
tions. The measured cross sections are generally found 
to agree with the NLO calculation although certain re- 
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FIG. 3: The ratio of pQCD predictions to the measured differ- 
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tainties due to variations of the factorization and renormal- 
ization scale. 



gions of phase space are identified where these predictions 
could better match the data. 

Supplementary material including tabulated W+n jet 
cross section measurements, theoretical predictions, and 
hadronization corrections applied to the theory can be 
found in the appendix to this document and online at 
doi:10.1016/j.physletb.2011. 10.011. 
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Appendix: Tables of Measurements, pQCD calculations and non-perturbative 
hadronization corrections 



In this appendix, we provide tables of the measured differential cross sections, theory predictions and hadronization 
corrections described in this paper. The region that defines the kinematic phase space of the measurement at particle 
level is given by the electron transverse momentum, pj, > 15 GeV, and pseudorapidity |?7 e | < 1.1, total transverse 
energy of all neutrinos Ej, > 20 GeV, W transverse mass Mjf > 40 GeV, jet transverse momentum p 3 ^ > 20 GeV 
and rapidity \y }ct \ < 3.2. W boson inclusive cross sections per jet multiplicity bin correspond to the sum over all pi!^ 
in the given jet multiplicity, and the normalized cross sections are the absolute cross sections in a given jet multiplicity 
divided by the inclusive W boson cross sections in the entire kinematic region. The hadronization corrections can 
be applied as a multiplicative factor to parton level jets clustered using the SISCONE algorithm to produce particle 
jets, as defined by the DO midpoint algorithm. Quoted systematic uncertainties on measured absolute cross-sections 
include a 6.1% luminosity uncertainty. 



TABLE I: Measured differential cross section, normalized to the measured inclusive W cross section, as a function of leading 
jet pr for events with one or more identified jets, along with statistical and systematic uncertainties. 

^-iLeading jet 



(GeV) 



^Leading jet ^ 



20-30 
30-40 
40-50 
50-60 
60-70 
70-80 
80-90 
90 - 110 
110 - 130 
130 - 150 
150 - 170 
170 - 210 
210 - 300 



24.6 

34.7 

44.7 

54.8 

64.8 

74.8 

84.8 

99.3 

119.3 

139.3 

159.2 

187.4 

241.8 



(5.95 ± 0.02(stat)l°;g2(syst) 
(2.30 ± 0.01(stat)l ( ( j'^(syst) 
(l.l3 ± O.Ol(stat) ± 0.04(syst) 
f 6.16 ± 0.05(stat) ± 0.26(syst) 

(3.43 ± 0.03(stat)tg;iij(syst) 
(2.OO ± 0.02(stat)tJ5;[J!(syst) 
M.22 ± 0.02(stat)lo;°s ( s y st ) 
(6.43 ± 0.09(stat) ±o.l(syst) 
(2.76 ± 0.05(stat)lg;^(syst) 
(l.29 ± 0.04(stat)lg;Ji(syst) 
(6.45 ± 0.24(stat)^'g(syst) 
(2.53 ± 0.10(stat)l C o|g(syst) 
( 4.48 ± 0.28(stat)±°'®°(syst) 



(Normalized) cross section: 



0.109 ± 0.0002(stat)to;^5(syst) 
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TABLE II: Normalized NLO theory predictions for W+l jet cross sections before application of hadronization corrections for 
the leading jet pr- 

g Loadin S jct l / aw . da / dp^ di " s iat (1/GcV) +1 standard deviation -1 standard deviation 



MCFM cross section predictions normalized by the MCFM inclusive W boson cross section 



20-30 


5.91 x 10 -3 


+3.0 x 10" 4 


-3.7 x 10~ 4 


30-40 


2.74 x 10~ 3 


+1.9 x 10~ 4 


-1.8 x 10~ 4 


40-50 


1.38 x 10~ 3 


+0.9 x 10~ 4 


-1.1 x 10~ 4 


50-60 


7.30 x 10~ 4 


+5.5 x 10~ 5 


-5.4 x 10~ 5 


60-70 


4.16 x 10~ 4 


+3.5 x 10~ 5 


-3.7 x 10~ 5 


70-80 


2.44 x 10~ 4 


+1.8 x 10~ 5 


-2.4 x 10~ 5 


80-90 


1.44 x 10~ 4 


+1.7 x 10~ 5 


-0.8 x 10~ 5 


90-110 


7.59 x 10~ 5 


+6.6 x 10~ 6 


-6.9 x 10~ 6 


110-130 


3.35 x 10~ 5 


+3.2 x 10 -6 


-3.5 x 10~ 6 


130-150 


1.56 x 10~ 5 


+1.1 x 10~ 6 


-1.6 x 10~ 6 


150-170 


7.34 x 10~ 6 


+1.3 x 10 -6 


-0.6 x 10~ 6 


170-210 


3.28 x 10~ 6 


+2.4 x 10~ 7 


-6.6 x 10~ 7 


210-300 


5.30 x 10~ 7 


+1.9 x 10~ 8 


-8.0 x 10" 8 


blackhat+SHERPA cross section predictions normalized by the BLACKHAT+SHERPA inclusive W boson cross section 


20-30 


6.14 x 10~ 3 


+1.9 x 10~ 4 


-3.1 x 10 -4 


30-40 


2.71 x 10~ 3 


+1.0 x 10~ 4 


-1.5 x 10~ 4 


40-50 


1.34 x 10~ 3 


+4.7 x 10~ 5 


-7.7 x 10~ 5 


50-60 


7.18 x 10~ 4 


+2.9 x 10~ 5 


-4.5 x 10~ 5 


60-70 


4.05 x 10~ 4 


+2.1 x 10~ 5 


-2.8 x 10~ 5 


70-80 


2.37 x 10~ 4 


+1.3 x 10~ 5 


-1.7 x 10~ 5 


80-90 


1.43 x 10~ 4 


+0.8 x 10" 5 


-1.0 x 10~ 5 


90-110 


7.49 x 10~ 5 


+5.3 x 10~ 6 


-6.0 x 10~ 6 


110-130 


3.31 x 10~ 5 


+2.6 x 10~ 6 


-2.9 x 10~ 6 


130-150 


1.55 x 10~ 5 


+1.3 x 10~ 6 


-1.4 x 10~ 6 


150-170 


7.43 x 10~ 6 


+6.0 x 10" 7 


-6.4 x 10~ 7 


170-210 


3.17 x 10~ 6 


+3.4 x 10~ 7 


-3.3 x 10~ 7 


210-300 


5.46 x 10~ 7 


+5.6 x 10~ 8 


-6.0 x 10~ 8 



TABLE III: Hadronization corrections derived with SHERPA 1.2.3 and the CTEQ6.6 PDF set for the leading jet pr distribution. 



4 c * dms]t ' (GeV) 


Hadronization correction 


20-30 


1.110 ± 0.004 


30-40 


0.968 ± 0.005 


40-50 


0.907 ± 0.005 


50-60 


0.914 ± 0.007 


60-70 


0.892 ± 0.008 


70-80 


0.898 ± 0.009 


80-90 


0.87 ± 0.01 


90-110 


0.87 ± 0.01 


110-130 


0.85 ± 0.01 


130-150 


0.85 ± 0.02 


150-170 


0.82 ± 0.02 


170-210 


0.79 ± 0.02 


210-300 


0.74 ± 0.03 
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TABLE IV: Measured differential cross section, normalized to the measured inclusive W cross section, as a function of second 
jet pr for events with two or more identified jets, along with statistical and systematic uncertainties. 



^Second jot CQpy| 


/^Second jotv fQ p y\ 


1 law ■ rf(r/f/n Socond jot (1 /GeVI 


20-30 


24.4 


('l.l0±0.01(stat)l ( ( J ,;^(syst)) x 10" 3 


30-40 


34.6 


(3.38±0.03(stat)l°;^(syst)J x 10~ 4 


40-50 


44.7 


(l.39±0.02(stat) + 2™(syst)) x 10~ 4 


50-60 


54.7 


(6.53 ± 0.12(stat) ± 0.37(syst)) x 10~ 5 


60-70 


64.7 


(3.33±0.09(stat)tg;i|(syst)J x 10~ 5 


70-80 


74.8 


M.80±0.06(stat)l°^(syst)J x 10~ 5 


80 - 100 


89.0 


(8.51 ±0.30(stat)lg;|o(syBt)J x 10~ 6 


100 - 140 


117.1 


(2.38±0.11(stat)l C o^(syst)J x 10~ 6 


140 - 180 


155.8 


(5.15±0.51(stat)l:|3; 4 ?(syst)j x 10~ 7 


180 - 250 


199.3 


(9.40 + 1.6(stat)l}; 3 (syst)J x 10~ 8 


(Normalized) cross section: 




0.017 ± 0.0001 (stat) +*-°* ^(syst) 



TABLE V: Normalized NLO theory predictions for W+2 jet cross sections before application of hadronization corrections for 
the second jet pr- 

gf, ccondjct (GeV) l/a w ■ da/dp^ cand ict (1/GeV) +1 standard deviation -1 standard deviation 



mcfm cross section predictions normalized by the mcfm inclusive W boson cross section 



20-30 


1.1 x 10" 3 


+6.2 x 10" 5 


-7.3 x 10~ b 


30-40 


3.6 x 10~ 4 


+5.1 x 10~ 5 


-1.5 x 10~ 5 


40-50 


1.5 x 10~ 4 


+0.0 


-1.4 x 10~ 5 


50-60 


7.1 x 10~ 5 


+4.9 x 10~ 6 


-5.7 x 10~ 6 


60-70 


3.3 x 10~ 5 


+0.0 


-2.8 x 10~ 6 


70-80 


2.0 x 10~ 5 


+5.4 x 10~ 6 


-1.7 x 10~ 6 


80-100 


8.1 x 10~ 6 


+1.0 x 10~ 6 


-8.1 x 10~ 7 


100-140 


2.7 x 10~ 6 


+0.0 


-2.9 x 10~ 7 


140-180 


4.3 x 10~ 7 


+2.4 x 10~ 7 


-3.9 x 10~ 8 


180-250 


7.5 x 10~ 8 


+0.0 


-3.4 x 10~ 9 


blackhat+SHERPA cross section predictions normalized by the blackhat+SHERPA inclusive W boson cross section 


20-30 


1.2 x 10~ 3 


+0.7 x 10~ 4 


-1.3 x 10" 4 


30-40 


3.7 x 10^ 4 


+2.1 x 10~ 5 


-4.1 x 10~ 5 


40-50 


1.5 x 10~ 4 


+0.7 x 10~ 5 


-1.5 x 10~ 5 


50-60 


7.1 x 10~ 5 


+3.9 x 10~ 6 


-7.6 x 10~ 6 


60-70 


3.6 x 10~ 5 


+1.2 x 10~ 6 


-3.5 x 10~ 6 


70-80 


2.0 x 10~ 5 


+0.8 x 10~ 6 


-1.8 x 10~ 6 


80-100 


9.4 x 10~ 6 


+3.5 x 10~ 7 


-9.1 x 10~ 7 


100-140 


2.6 x 10~ 6 


+0.4 x 10~ 7 


-2.2 x 10~ 7 


140-180 


5.5 x 10~ 7 


+0.0 


-4.6 x 10~ 8 


180-250 


1.0 x 10~ 7 


+0.0 


-1.0 x 10~ 8 
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TABLE VI: Hadronization corrections derived with sherpa 1.2.3 and CTEQ6.6 PDF set for the second jet pr distribution. 



g sccondjct ( GeV ) Hadronization correction 

20-30 1.133 ± 0.006 

30-40 0.942 ± 0.009 

40-50 0.89 ± 0.01 

50-60 0.92 ± 0.02 

60-70 0.88 ± 0.02 

70-80 0.82 ± 0.03 

80-100 0.85 ± 0.03 

100-140 0.88 ± 0.03 

140-180 0.85 ± 0.05 

180-250 0.62 ± 0.05 



TABLE VII: Measured differential cross section, normalized to the measured inclusive W cross section, as a function of third 
jet pr for events with three or more identified jets, along with statistical and systematic uncertainties. 



£Thirdjet (GeV) 


(£™ rdicl ) (GeV) 


l/a w ■ da/dp™' d Jot (1/GeV) 


20-30 


24.3 ( 


'2.00±0.04(stat)lg;l|(syBt)) x 10" 4 


30-50 


38.1 ( 


2.89±0.10(stat)t°'j«(syst)J x 10~ 5 


50-80 


61.8 ( 


2.54±0.20(stat)t°'j 4 (syst)J x 10~ 6 


80 - 130 


96.2 | 


1.48±0.26(stat)±°'j 4 (syst)J x 10~ 7 



(Normalized) cross section: 0.0026 ± 0.5 x 10" 4 (stat)+|^^_* (syst) 



TABLE VIII: Normalized NLO theory predictions before application of hadronization corrections for third jet pr- 

ff™ rdict (GeV) l/aw ■ da/dpl h " dict (1/GeV) +1 standard deviation -1 standard deviation 

ROCKET+MCFM cross section predictions normalized by the rocket+mcfm inclusive W boson cross section 



20-30 


1.7 x 10~ 4 


+5.3 x W~ b 


-3.1 x 10~ b 


30-50 


2.7 x 10~ 5 


+4.2 x 10~ 6 


-5.2 x 10~ 6 


50-80 


2.2 x 10~ 6 


+3.0 x 10~ 7 


-3.6 x 10~ 7 


80-130 


1.2 x 10~ 7 


+0.7 x 10" 8 


-2.4 x 10~ 8 



BLACKHAT+ sherpa cross section predictions normalized by the blackhat+SHERPA inclusive W boson cross section 
20-30 1.8 x 10~ 4 +2.1 x 10~ s -3.2 x 10" a 

30-50 2.6 x 10~ 5 +3.0 x 10~ 6 -4.3 x 10~ 6 

50-80 2.2 x 10~ 6 +2.4 x 10~ 7 -3.6 x 10~ 7 

80-130 1.3 x 10~ 7 +1.2 x 10~ 8 -2.1 x 10~ 8 



TABLE IX: Hadronization corrections derived with sherpa 1.2.3 and the CTEQ6.6 PDF set for the third jet pr distribution. 



£Thirdjet ^ QeY) 


Hadronization correction 


20-30 


1.09 ± 0.01 


30-50 


0.87 ± 0.01 


50-80 


0.83 ± 0.02 


80-130 


0.75 ± 0.03 
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TABLE X: Measured differential cross section, normalized to the measured inclusive W cross section, as a function of fourth 
jet pr for events with four or more identified jets, along with statistical and systematic uncertainties. 



^Fourth jet ( Ge y) 


^Fourth jet ^ (GeV) 


1/W ■ da/dp^° urth Jct (1/GeV) 


20-30 


24.1 


(3.00±0.17(stat)i°;^(syst)') x 10~ 5 


30-50 


37.6 


(2.46±0.31(stat)l°-j?(syst)J x 10~ 6 


50 - 80 


60.3 


(9.87±2.78(stat)t|;^(syst)J x 10~ 8 


(Normalized) cross section: 




0.00035 ± 1.8 x 10" 5 (stat)+^ '^^^(syst) 



TABLE XI: Normalized LO theory predictions before application of hadronization corrections for fourth jet pr- 
g Fo„ rth g ( GeV ) l j (Tw ; rft7 / dp F°urth jet (j/Q^) +1 standard deviation -1 standard deviation 

ROCKET+MCFM cross section predictions normalized by the rocket+mcfm inclusive W boson cross section 
20-30 2.9 x 10 -5 +2.5 x 10~ b -1.2 x 10 _b 

30-50 2.7 x 10~ 6 +2.3 x 10~ 6 -1.1 x 10~ 6 

50-80 1.2 x 10~ 7 +1.1 x 10~ 7 -5.4 x 10~ 8 

blackhat+SHERPA cross section predictions normalized by the blackhat+SHERPA inclusive W boson cross section 
20-30 3.0 x 10 _b +2.5 x 10~ b -1.2 x 10 _b 

30-50 2.7 x 10~ 6 +2.3 x 10~ 6 -1.1 x 10~ 6 

50-80 1.2 x 10~ 7 +9.8 x 10~ 8 -4.7 x 10~ 8 



TABLE XII: Hadronization corrections derived with SHERPA 1.2.3 and the CTEQ6.6 PDF set for the fourth jet pr distribution. 



£ F ourthj 


et (GcV) 


Hadronization correction 


20-30 




1.12 ± 0.02 


30-50 




0.86 ± 0.02 


50-80 




0.74 ± 0.03 
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TABLE XIII: Measurement of n th jet inclusive cross section in association with a 


W boson. All results are given in pb. 


Bin 


a w+n-\ct (pb) Statistical Uncertainty (pb) 


Systematic Uncertainty (pb) 


W + 0-jet 
W + 1-jet 
W + 2-jet 
W + 3-jet 
W + 4-jet 


1097 ± 0.8 (stat) 
119.5 ± 0.3 (stat) 
19.0 ± 0.1 (stat) 
2.9 ±0.1 (stat) 
0.39 ± 0.02 (stat) 


-89 i b y bi, j 
ts.s (syst) 
t?:3 (^st) 
±8:1 (syst) 
i°o:°o? (syst) 


TABLE XIV: blackhat+SHERPA and ROCKET+MCFM predictions of n th jet inclusive cross sections in association with a W 
boson. All results are given before hadronization corrections. 


Bin 


CT w+n-jet (P b ) 


Scale Uncertainty (pb) 


BLACKHAT+SHERPA, no hadronization correction 


W + 0-jet 
W + 1-jet 
W + 2-jet 
W + 3-jet 
W + 4-jet 


1153 
138.1 
21.4 
2.81 
0.31 


t? (syst) 
±l\ (syst) 
±1:1 (syst) 
±o.5o (syst) 
i'o 1 :^ (syst) 


rocket+mcfm, no hadronization correction 


W + 0-jet 
W + 1-jet 
W + 2-jet 
W + 3-jet 
W + 4-jet 


1207 
143.6 
21.6 
2.82 
0.32 


1? (syst) 

tio.l (syst) 
±11 (syst) 

±0.53 (Syst) 

±S:iI (syst) 
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TABLE XV: Measurement of n th jet to (n — l) th jet inclusive cross section ratios, 
ffn/fn-i Ratio (data) 

a 2 /a, 0.159±°;8g| 



0-3/0-2 0.154±°- 
0-4/0-3 0. 132±g;gii 



005 
006 



TABLE XVI: blackhat+SHERPA and rocket+mcfm predictions of n th jet to (n - l) th jet inclusive cross section ratios before 
hadronization corrections. 



<7 n /(Tn-l 


Ratio (blackhat+sherpa, 


no hadronization correction) 


0"1 1 (To 
O-2/0-I 
O-3/C2 
O-4/0-3 




120 +UUU4 

U.lziU_ 007 
n lcc+0.003 
U.±OO_ oog 

n 101 +0.007 
u,ioi -0 .009 

u - 11 -0.01 


o n ja n -\ 


Ratio (rocket+mcfm, 


no hadronization correction) 


0-1/0-0 
0-2/0-1 
0-3/0-2 
0-4/0-3 




q n g+0.008 
U.ll»_ 008 
n 1 ci +0.0005 
U.101_o.oil 

U.lOU_ 007 

U.10_o.oi 



TABLE XVII: Inclusive hadronization correction derived with SHERPA 1.2.3 and the CTEQ6.6 PDF set for each inclusive jet 
bin. 



Inclusive jet multiplicity bin Hadronization correction and associated statistical uncertainty 

0- jet LOO" 

1- jet 1.021 ± 0.003 

2- jet 1.052 ± 0.006 

3- jet 1.035 ± 0.011 

4- jet 1.078 ± 0.022 



